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Abstract: We consider leptogenesis induced by soft supersymmetry breaking terms ("soft 
leptogenesis"), in the context of the inverse seesaw mechanism. In this model there are 
^ I lepton number (L) conserving and L-violating soft super symmetry-breaking B-terms in- 

volving the singlet sneutrinos which, together with the - generically small- L-violating 
parameter responsible of the neutrino mass, give a small mass splitting between the four 
singlet sneutrino states of a single generation. In combination with the trilinear soft super- 
symmetry breaking terms they also provide new CP violating phases needed to generate a 
lepton asymmetry in the singlet sneutrino decays. We obtain that in this scenario the lep- 
ton asymmetry is proportional to the L-conserving soft supersymmetry-breaking i?-term, 
and it is not suppressed by the L-violating parameters. Consequently we find that, as in the 
standard see-saw case, this mechanism can lead to sucessful leptogenesis only for relatively 
small value of the relevant soft bilinear coupling. The right-handed neutrino masses can be 
sufficiently low to elude the gravitino problem. Also the corresponding Yukawa couplings 
involving the lightest of the right-handed neutrinos are constrained to be ^ l^ifcP ^ 10"'' 
which generically implies that the neutrino mass spectrum has to be strongly hierarchical. 
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1. Introduction 

The discovery of neutrino oscillations makes leptogenesis a very attractive solution to the 
baryon asymmetry problem [1]. In the standard framework it is usually assumed that the 
tiny neutrino masses are generated via the (type I) seesaw mechanism [2] and thus the new 
singict neutral leptons with heavy (lepton number violating) Majorana masses can produce 
dynamically a lepton asymmetry through out of equilibrium decay. Eventually, this lepton 
asymmetry is partially converted into a baryon asymmetry due to fast B — L violating 
sphaleron processes. 

For a hierarchical spectrum of right-handed neutrinos, successful leptogenesis requires 
generically quite heavy singlet neutrino masses [3], of order M > 2.4(0.4) x 10^ GeV 
for vanishing (thermal) initial neutrino densities [3,4], although flavour effects [5] and/or 
extended scenarios [6] may affect this limit *. The stability of the hierarchy between this 
new scale and the electroweak one is natural in low-energy supersymmetry, but in the 
supersymmetric seesaw scenario there is some conflict between the gravitino bound on the 
reheat temperature and the thermal production of right-handed neutrinos [9]. This is so 
because in a high temperature plasma, gravitinos are copiously produced, and their late 
decay could modify the light nuclei abundances, contrary to observation. This sets an 
upper bound on the reheat temperature after inflation, Trh < 10^"^^ GeV, which may be 
too low for the right-handed neutrinos to be thermally produced. 

Once supersymmetry has been introduced, leptogenesis is induced also in singlet sneu- 
trino decays. If supersymmetry is not broken, the order of magnitude of the asymmetry and 

*This bound applies when the lepton asymmetry is generated in the decay of the lightest right-handed 

neutrino. The possibility to evade the bound producing the asymmetry from the second lightest right- 
handed neutrino has been considered in [7], and flavour effects have been analysed for this case in [8]. 
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the basic mechanism are the same as in the non-supersymmetric case. However, as shown 
in Refs. [10, 11], supersymmetry-breaking terms can play an important role in the lepton 
asymmetry generated in sneutrino decays because they induce effects which are essentially 
different from the neutrino ones. In brief, soft supersymmetry-breaking terms involving 
the singlet sneutrinos remove the mass degeneracy between the two real sneutrino states of 
a single neutrino generation, and provide new sources of Icpton number and CP violation. 
As a consequence, the mixing between the two sneutrino states generates a CP asymme- 
try in the decay, which can be sizable for a certain range of parameters. In particular 
the asymmetry is large for a right-handed neutrino mass scale relatively low, in the range 
10^ — 10^ GeV, well below the reheat temperature limits, what solves the cosmological 
gravitino problem. Moreover, contrary to the traditional leptogenesis scenario, where at 
least two generations of right-handed neutrinos are required to generate a CP asymmetry 
in neutrino/sneutrino decays, in this new mechanism for leptogenesis the CP asymmetry 
in sneutrino decays is present even if a single generation is considered. This scenario has 
been termed "soft leptogenesis", since the soft terms and not flavour physics provide the 
necessary mass splitting and CP-violating phase. It has also been studied in the minimal 
super symmetric triplet seesaw model [12]. 

In this paper we want to explore soft leptogenesis in the framework of an alternative 
mechanism to generate small neutrino masses, namely the inverse seesaw scheme [13]. This 
scheme is characterized by a small lepton number violating Majorana mass term fi, while 
the effective light neutrino mass is rrii^ oc fi. Small values of ji are technically natural, 
given that when /x — ^ a larger symmetry is realized [14] : lepton number is conserved and 
neutrinos become massless. In the inverse seesaw scheme lepton flavour and CP violation 
can arise even in the limit where lepton number is strictly conserved and the light neutrinos 
are massless [15], due to the mixing of the SU(2) doublet neutrinos with new SU(2) x U(l) 
singlet leptons. 

As opposite to the standard seesaw case, these singlet leptons do not need to be 
very heavy [16], and, as a result, lepton flavour and CP violating processes are highly 
enhanced [15]. In Ref. [17] it was studied the possibility that the baryon asymmetry is 
generated in this type of models during the electroweak phase transition, in the limit 
/Lt = 0. A suppression was found due to the experimental constraints on the mixing angles 
of the neutrinos [18]. Therefore we consider here the supersymmetric version of the model 
and the soft leptogenesis mechanism, since {i) in this case we expect that a CP asymmetry 
will be generated in sneutrino decays even with a single-generation and no suppression due 
to the mixing angles is expected, and (ii) this scheme provides a more natural framework 
for the relatively low right-handed neutrino mass scale. 

The outline of the paper is as follows. Section 2 presents the main features of the 
inverse seesaw model in the presence of supersymmetry breaking terms. In Sec. 3 we 
evaluate the lepton asymmetry generated in the decay of the singlet sneutrinos using a 
field-theoretical approach assuming a hierarchy between the SUSY and L breaking scales 
and the mass scale of the singlet sneutrinos. The relevant Boltzmann equations describing 
the decay, inverse decay and scattering processes involving the singlet sneutrino states are 
derived in Sec. 4. Finally in Sec. 5 we present our quantitative results. In appendix A 
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we recompute the asymmetry using a quantum mechanics approach, based on an effective 
(non hermitic) Hamiltonian [10,11]. 



2. Inverse Seesaw Mechanism 

In this type of models [13], the lepton sector of the Standard Model is extended with two 
electroweak singlet two-component leptons per generation, i.e., 

L'=(^2y^R^^R^''i^ (2.1) 

We assign lepton number L = 1 to the singlets s^^ and Uj^. In the original formulation 
of the model, the singlets were superstring inspired E(6) singlets, in contrast to the 
right-handed neutrinos Vj^, which are in the spinorial representation. More recently this 
mechanism has also arisen in the context of left-right symmetry [19] and SO(IO) unified 
models [20]. 

The (9 X 9) mass matrix of the neutral lepton sector in the ul, f^, sl basis is given by 



M = 



/ mn 

I (2.2) 

^ M fi 



where mo^M are arbitrary 3x3 complex matrices in flavour space and /i is complex sym- 
metric. In models where lepton number is spontaneously broken by a vacuum expectation 
value (fj), /X = \{a) [21]. The matrix M. can be diagonalized by a unitary transforma- 
tion, leading to nine mass eigenstates Ua- three of them correspond to the observed light 
neutrinos, while the other three pairs of two component leptons combine to form three 
quasi-Dirac leptons. 

In this "inverse seesaw" scheme, assuming m^),^ <C M the effective Majorana mass 
matrix for the light neutrinos is approximately given by 

= m'^M'^~^ iJ.M~^mD , (2.3) 

while the three pairs of heavy neutrinos have masses of order M, and the admixture among 
singlet and doublet SU{2) states is suppressed by mo/M. Although M is a large mass 
scale suppressing the light neutrino masses, in contrast to the Majorana mass (AL = 2) of 
the right-handed neutrinos in the standard seesaw mechanism, it is a Dirac mass (AL = 0), 
and it can be much smaller, since the suppression in Eq. (2.3) is quadratic and moreover 
light neutrino masses are further suppressed by the small parameter jj, which characterizes 
the lepton number violation scale. 

Notice that in the ^ limit lepton number conservation is restored. Then, the 
three light neutrinos are massless Weyl particles and the six heavy neutral leptons combine 
exactly into three Dirac fcrmions. 

We arc going to consider the supcrsymmetric version of this model (for an analysis of 
lepton flavour violation in this case see [22]). In this case, the above neutral lepton mass 
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matrix (2.2) is described by the following superpotential: 




3 ' 



(2.4) 



where i^j = 1,2,3 are flavour indices, i?, Lj, A^j, Si are the superfields corresponding to the 
SU (2) up-Higgs and lepton doublets, and i^)^ and s\ singlets, respectively, and Yij denote 
the neutrino Yukawa couplings. Thus, after spontaneous electroweak symmetry breaking, 
the neutrino Dirac masses are given by 



The relevant soft supcrsymmctry breaking terms arc the bilinear and trilinear scalar cou- 
plings involving the singlet sneutrino fields, that provide new sources of lepton number 
and CP violation. Prom now on we consider a simplified one generation of {N, S) model 
because a single generation of singlet sneutrinos is sufficient to generate the CP asymme- 
try. Indeed in the three-generation case, the relevant out of equilibrium decays are usually 
those of the lightest heavy singlet states while the decay of the heavier (if heavier enough) 
give no effect. Thus with our simplified single generation model we refer to the lightest of 
the three heavy singlet sneutrinos which we number as 1. Consequently we label M = Mn 
and )U = Also, for simplicity, we will assume proportionality of the soft trilinear terms. 



-Lsoft = AYuLiNH + ih^SS^ + m%NN^ + m|^5iV^ + BsSS + BsnSN + h.c. (2.6) 



With our lepton number assignments, the soft SUSY breaking terms which violate L are 
m|jY and Bs- The sneutrino interaction Lagrangian is then: 



This Lagrangian has three independent physical CP violating phases: (f)B which can be 
assigned to Bsn, <i>A which is common to the three terms with AYn^ and (f)MY which is 
common to theT three terms with MY^^, and are given by: 



where we have defined = ^uM* + m'gj^. These phases provide the CP violation 

necessary to generate dynamically a lepton asymmetry, even with a single generation of 
sneutrinos. They can also contribute to lepton electric dipole moments [23]. 



{mD)ij = Yij{H) 



(2.5) 



L = {YuLiNH + YuLiNh + YuLiNh + h.c.) + 
+ iYuM*LiS^H + AYuLiNH + h.c.) + 
+ {fiSS + MSN + h.c.) + 

+ ((1^1^ + m| + \Mf)S^S + {m% + \M\^)NN^ + {fxM* + m|jv)5iV^ + h.c.) + 
+ {BsSS + BsnSN + h.c.) 



(2.7) 



cPb = wg{BsNB*sMj^) 
<Pa = arg(AB* MV*(M|^)2) 
(f)MY = arg(M|;|VM» , 



(2.8) 
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(2.9) 



From the Lagrangian in Eq. (2.7) we obtain the sneutrino mass matrix in the interaction 
basis {Fi = N,N^,S,S^): 

fm% + \M\^ M|*^ BsN \ 

m% + \M\' B*s^ Ml^ 

MIj^ Bsn |//|' + m| + |M|2 2Bs 

V Mf^ 2B*s \^lf + ml + \M\'J 

Notice that in the most general case it is not possible to remove all CP phases from the 
sneutrino mass matrix. With our choice of basis, Bsn is the only complex parameter, 
Bsn = \BsN\e''f"' . 

Although one can easily obtain the analytic expressions for the corresponding mass 
eigenvalues and eigenvectors, for the general case they arc lengthy and we do not give 
them here. Under the assumption that all the entries are real, i.e., = 0, we obtain the 
following mass eigenvalues: 

Mf = M'2 + 55 - ^ ^J4{BsN + Ml^f + (2^5 - + m| + \^Ji\^f 

1 
2 



M| = M'2 -Bs + -\l^{BsN - M|^)2 + (255 + - m| - |//|2)2 



M| = M'2 -Bs- \\J^S!s - Mh? + (2i?,s' + ml - - 



Ml = M'^ + Bs + \ ^l^{BsN + M|^)2 + (2^5 -ml+fhl + , (2.10) 

where we have defined M'^ = |Mp + m% + ml + 

Furthermore, if we assume conservative values of the soft breaking terms: 



A ~ 0{msusY) 
ffiN ^ rhs ^ fhsN ~ 0{msusY) 
Bs ~ 0{msuSYli) 
Bsn ~ 0{msuSYM) 



(2.11) 



with both, n,msuSY ^ M, we see that Bs,m1f,'m'g,m'g^ <C Bsn and Mg^ ~ /xM*. 
Neglecting these small soft terms, there is still one physical CP violating phase, 

<i) = (l)A-^B = argiAB*sNM) . (2.12) 

In this limit, we choose for simplicity a basis where A = l^je*"^ is the only complex param- 
eter. Then we diagonalize to first order in two expansion parameters, 

ImI ~ \Bsn\ 



2\M\ ' 2|Mp 
To this order the mass eigenvalues are: 



0{msuSY/M) (2.13) 



Mf = M"^ - Mil- BsN 
M| = M"^ - Mfi + BsN 

M| = M^ + Mfi- BsN (2.14) 
M| = M'^ + Mii + Bsn 
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and the eigenvectors: 



Ni = 


l(s.-St) 




N2 = 






N3 = 






7V4 = 


1 + 





(2.15) 



Note that in this limit the mass degeneracy among the four sneutrino states is removed 
by both the L-violating mass fi and L-conserving supersymmctry breaking term Bsn- 
Together with the trihnear A term they also provide a source of CP violation, and the 
mixing among the four sneutrino states leads to a CP asymmetry in their decay. 

Another interesting limit is to diagonalize the sneutrino mass matrix (2.9) neglecting 
only the Bs entry, which may be appropriate if ^ <C msuSY and the order of magnitude 
of the soft breaking terms is as given by (2.11). In this limit the mass matrix can also 
be taken real and the mass eigenvalues can be read from Eq. (2.10), just setting Bs = 0. 
Now there are two non zero CP violating phases, (pYM and 0'^ = aj:g{ABgj^M'^iJ,*Mgj^). 
However the combination that is relevant for the CP asymmetry in sneutrino decays is the 
same as in the previous case, (f) = ^ym + (t>'A~ ^^si^Bg^M). 

As we will see in Sec. 5, the total CP asymmetry in the singlet sneutrino decays turns 
out to be sizable for very small values of the soft term Bsn ^ MmsusY- Neglecting the 
Bsn term in the Lagrangian there are still two CP violating phases, (pMY and but 
again the sneutrino mass matrix can be taken real, so that the mass eigenvalues are as 
given by Eq.(2.10) with Bsn = 0. The phase relevant for the CP asymmetry in the singlet 
sneutrino decays is now (f)' = (f)YM + ^A = &^g{AB*gM Mg^) . 

Finally, neglecting supersymmetry breaking effects, the total singlet sneutrino decay 
width is given by 



3. The CP Asymmetry 

In this section we compute the CP asymmetry in the singlet sneutrino decays. As discussed 
in Ref. [11], when V S> AMjj = Mj — Mj, the four singlet sneutrino states are not well- 
separated particles. In this case, the result for the asymmetry depends on how the initial 
state is prepared. In what follows we will assume that the sneutrinos are in a thermal bath 
with a thermalization time r"-"^ shorter than the typical oscillation times, AM^^, therefore 
coherence is lost and it is appropriate to compute the CP asymmetry in terms of the mass 
eigenstates Eq.(2.15). 
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The CP asymmetry produced in the decay of the state iVj is given by (see sec. 4): 

^r(iv,->/)-r(iv,^/) 

= 4^ ^ ^ , (3.1) 

f 

where / = L^H, Lkh. We also define the fermionic and scalar CP asymmetries in the decay 
of each iVj as 



_k 

k 

^ \Ai{Ni ^ Lfc/i)|2 - \Ai{Ni ^ Zfe/i)|2 



(3.2) 



Notice that Cg. and e/. are defined in terms of decay amphtudes, without the phase-space 
factors which, as we will see, are crucial to obtain a non-vanishing CP asymmetry, much 
as in the standard seesaw case [10, 11]. The total asymmetry Cj generated in the decay of 
the singlet sneutrino iVj can then be written as 

= ^-El±lf^ , (3.4) 

Cs + Cf 

where Cs,Cf are the phase-space factors of the scalar and fermionic channels, respectively. 

Since the scale of lepton number and supersymmetry breaking are fi, msuSY ^ M, 
there is an enhancement of the CP violation in mixing (wave-function diagrams) , so we only 
include this leading effect and neglect direct CP violation in the decay (vertex diagrams). 

We compute the CP asymmetry following the effective field theory approach described 
in [25], which takes into account the CP violation due to mixing of nearly degenerate 
states by using resumed propagators for unstable (mass eigenstate) particles. The decay 
amplitude of the unstable external state Ni defined in Eq. (2.15) into a final state / is 
described by a superposition of amplitudes with stable final states: 

AJN.^ f) = A{ -yA{ , '^l' , (3.5) 

where A{ are the tree level decay amplitudes and Iljj are the absorptive parts of the two- 
point functions for i,j = 1,2,3,4. The amplitude for the decay into the conjugate final 
state is obtained from (3.5) by the replacement Aj A^*. 
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The decay amplitudes can be read off from the interaction Lagrangian (2.7), after 
performing the change from the current to the mass eigenstate basis: 

+ '-NsiYikLkh - {YikM* - AYik)LkH] 

+ ^N^Y.kLkh + {Y,kM* + AY,k)LkH] + h.c. (3.6) 

Neglecting supersymmetry breaking in vertices, up to an overall normalization we 
obtain that the decay amplitudes into scalars, A^*" = A{Ni L^H), verify A2'' = A3'' = 
iA^'' = iA\. Correspondingly the decay amplitudes into fermions = A{Ni Lj^h), 
verify A^'' = -A^* = -iA('' = lA^^. 

Keeping only the lowest order contribution in the soft terms, 

Uii = MY i = 1,...,4 (3.7) 

ni2 = Hsi = -n34 = -n43 = |a| r sinc^ (3.8) 

Altogether we can then write the fermionic and scalar CP asymmetries as: 

2(M2 - M])Tlji Im(A''*^f ) 

y — — (3-9) 



k 

fk* /ifk^ 



2{Mf - Mf)Uji J2 Im(A/**A 
k 

^ [{Mi-Mff + U%]J2\A: 



k 

Inserting the values of the amplitudes A^ and the absorptive parts of the two-point func- 
tions (3.7) we obtain the final expression for the scalar and fermionic CP asymmetries at 
T = 0: 

4|I?^vv A|r . 
= -'^^ = '^ = - 4|^5ArP + |Mpr^ ""^ 

and the total CP asymmetry generated in the decay of the sneutrino Ni is then: 

6,(r) = e, . (3.12) 

As long as we neglect the zero temperature lepton and slepton masses and small Yukawa 
couplings, the phase-space factors of the final states are flavour independent. After includ- 
ing finite temperature effects they are given by: 

Cf = {l-XL- Xh)\{l,XL.Xh) [1 - /f] [1 - (3.13) 
Cs = A(l, XH, xj) [1 + fj^] [1 + ff\ (3.14) 
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where 

Ql = eME^\/T] - 1 

are the Bose-Einstein and Fermi-Dirac equihbrium distributions, respectively, and 

EL,h = ^:{l + XL,h-Xh,L), ^^^,Z = f (l + ^/f,L-^Z,J (3-17) 
A(l,x,y) = V(l + x-y)2-4x, x„ = (3.18) 

The thermal masses for the relevant supersymmetric degrees of freedom are [26]: 

ml{T) = 2ml{T) = gc/i + ^5?- + ^A?) , (3.19) 

m|(T) = 2mi(r) = (^,i + i,^)T^ (3.20) 

Here §2 and (^y are gauge couplings and At is the top Yukawa, renormalized at the appro- 
priate high-energy scale. 

As we will see in the next section, from Eq. (4.46), if the initial distributions of all four 
states iVj are equal, their total contribution to the total lepton number can be factorized 

as: _ _ _ 

e(T) = y -4 . (3.21) 

f 

Several comments are in order. We find that this leptogenesis scenario presents many 
features analogous to soft leptogenesis in seesaw models [10-12]: (i) The CP asymmetry 
(3.12) vanishes if Cg = cj, because then there is an exact cancellation between the asymme- 
try in the fermionic and bosonic channels. Finite temperature effects break supersymmetry 
and make the fermion and boson phase-spaces different Cg 7^ c/, mainly because of the fi- 
nal state Fermi blocking and Bose stimulation factors, {ii) It also displays a resonance 
behaviour: the maximum value of the asymmetry is obtained for 2Bsn/M ~ F. (iii) 
The CP asymmetry is due to the presence of supersymmetry breaking and irremovable CP 
violating phases, thus it is proportional to |i?5jvy4| sincj). 

As seen from Eq.(3.11) we obtain that the CP asymmetry is not suppressed by the 
lepton number violating scale fx. This may seem counterintuitive. However if /x = the 
four sneutrino states are pair degenerate, and we can choose a lepton number conserving 
mass basis, made of the {L = 1) states 



N^ = ±(^S^ + n) (3.22) 
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and their hermitian conjugates, with L = — 1, N'^ ■ Although there is a CP asymmetry 
in the decay of these sneutrinos, it is not a lepton number asymmetry (since in the hmit 
H = total lepton number is conserved) but just a redistribution of the lepton number 
stored in heavy sneutrinos and light lepton and slepton SU{2) doublets. At very low 
temperatures, T ^ M, when no heavy sneutrinos remain in the thermal bath, all lepton 
number is in the light species and obviously if we started in a symmetric Universe with no 
lepton number asymmetry it will not be generated. 

In other words, the total lepton number generated for the case with no lepton number 
violation is zero but it cannot be recovered by taking the limit /l( ^ of Eq.(3.11) because 
in the derivation of Eq.(3.11) it is assumed implicitly that the four singlet sneutrino states 
are non-degenerate and consequently it is only valid if /j. (or some of the other L violating 
parameters) is non zero. 

In appendix A we recompute the asymmetry using a quantum mechanics approach, 
based on an effective (non hermitic) Hamiltonian [10,11], and we get the same parametric 
dependence of the result, which differs only by numerical factors. Both expressions agree 
in the limit T < \Bsn/M\. 

As discussed at the end of the previous section there may be other interesting ranges 
of parameters beyond Eq.(2.13). Thus in order to verify the stability of the results to 
departures from this expansion wc have redone the computation of the CP asymmetry 
keeping all the entries in the sneutrino mass matrix, and just assuming that it is real. The 
expressions are too lengthy to be given here but let us simply mention that we have found 
that, in the general case, the CP asymmetries generated in the decay of each of the four 
singlet sneutrino states are not equal but they can always be written as: 

where the functions fi verify that for <C \Bsn\ and to any order in \Bs\: 

J2 fi{Bs, fi, MIj,) oc \Bsn\ . (3.24) 

i 

In the limiting case I-Bsat] <C I-B5I, m|^^y, |^^| the dominant term in the CP asym- 
metry at leading order in \Bs\ ~ \M\T <C is 

% I O IS I 

It also exhibits a resonant behaviour, described now by |55|r/(4|i?5p + jMpF^)^, how- 
ever the total CP asymmetry in this limit is further suppressed by a factor of order 

{\l^\mlusY)l\Mh\- 

Finally, let's comment that in the previous derivation we have neglected thermal cor- 
rections to the CP asymmetry from the loops, i.e., we have computed the imaginary part 
of the one-loop graphs using Cutkosky cutting rules at T = 0. These corrections are the 

same for scalar and fermion decay channels, since only bosonic loops contribute to the 
wave-function diagrams in both cases, so they are not expected to introduce significant 
changes to our results. 
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4. Boltzmann Equations 



We next write the relevant Boltzman equations describing the decay, inverse decay and 
scattering processes involving the sneutrino states. 

As mentioned above we assume that the sneutrinos are in a thermal bath with a ther- 
malization time shorter than the oscillation time. Under this assumption the initial states 
can be taken as being the mass eigenstates in Eq. (2.15) and we write the corresponding 
equations for those states and the scalar and fermion lepton numbers. The CP fermionic 
and scalar asymmetries for each TVj defined at T = are those given in Eq. (3.11). 

Let's notice that the CP asymmetries as defined in Eq. (3.11) verify e^- = — ej- = e,. 
However in order to better trace the evolution of the scalar and fermion lepton numbers 
separately we will keep them as two different quantities in writing the equations. 

Using CPT invariance and the above definitions for the CP asymmetries and including 
all the multiplicative factors we have: 

- LkH)? = ^\A{LiH^ ATOP ^ ^-^^\Ar\' , 

k k k 

^|i(iV, ^ Zt^t)|2 = Y^lMLkH - , 



k k k 



(4.1) 



where 



k k 

The Boltzmann equations describe the evolution of the number density of particles in the 
plasma: 

'^''^ + 2,Hnx = Yl ^fl:: [/^/- • • • (1 ± ± /i) • • • ^(^"^ ■■■^xj...y 



dt 

j,l,m 



where, 



^\m.. 



fxfj...{l±fi){l±fm)...W{X3---^lm...)] 



d^Px f d^pj f d^pi f d'^p. 



{27Tf2Ex J {2^f2Ej"'J {27rf2EiJ {2n)^2E,, 



and W{lm ■ ■ ■ ^ Xj . . . ) is the squared transition amplitude summed over initial and final 
spins. In what follows we will use the notation of Ref. [24] and we will assume that the 
Higgs and higgsino fields are in thermal equilibrium with distributions given in Eqs. (3.15) 
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and (3.16) respectively, while the leptons and sleptons are in kinetic equilibrium and we 
introduce a chemical potential for the leptons, and sleptons, /x^: 

1 



f-L 

h 
/zt 



1 

" exp[(EL-/xy)/r] + l ' 
1 

1 



(4.3) 



exp[(^j - i2s)/T] - 1 ■ 

Furthermore in order to eliminate the dependence in the expansion of the Universe we 
write the equations in terms of the abundances Yx, where Yx = rix/s. 

We are interested in the evolution of sneutrinos Y^, , and the fermionic Yc and scalar 
Y-^ lepton number, defined as Yc = (Yl — Yi)/2, Y-^ = (Y^ — Y-j^^)/2. The number density 
of sneutrinos is regulated through its decays and inverse decays, defined in the D — terms, 
while to compute the evolution of the fermionic and scalar lepton number we also need to 
consider the scatterings where leptons and sleptons are involved. The scattering terms are 
defined in the S — terms. 



dYi 



Ni 



dt 

dYc 
dt 

dY^ 
dt 



-Di -Di-Di-D, 



E (^^ - ^0 - 2^ - ^LZt + ^LZt - s^z + Sa 

i 

E (a - 51) -2S- S^~^, + S^~^, + S^~^ - S^^ 



(4.4) 
(4.5) 

(4.6) 



where 



k 



(4.7) 



k 



(4.8) 



k 



(4.9) 
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A; 



.5 = Ai 



sS = A. 



/z/^'(i + /zt)(i + ^ 



k,k' 



-/Zt/^'(1 + /Z)(l + rH')T.\Msub{LlH^ - L^'H)f 



k,k' 



12 

Lit - ^'■34 



sS'rTt = A. 



/L/r(l + /Zt)(l + fw)\Y.\Msub{Lkh ^ 4,i/t^ 



-/zt/^'(i - - /r)Ei^-f(4^^ - L,,h)\' 



k,k' 



k,k' 



12 



sSrT = A 



/L/r(l + /z)(l + /h')EI^-^'(^^^ ^ Lk>H)f- 



k,k' 



„c „ _ Al2 



/L/r(i + /zt)(i + /^')Ei^-''(^'^^ 

-/zt/^'(i - - /DEl^-^'l^I^^ ^ 



The S — terms are defined in terms of subtracted amplitudes, since the on-shell contribution 
is already taken into account through the decays and inverse decays in the D — terms. So 
for example: 



(4.17) 



where, 



\Mos{Lkh ^ Lk'h)\^ = A{Lkh^N, 



2 7r5{s — ITlj;^ 



(4.18) 



In writing Eqs.(4.4)-(4.6) we have not included the AL = 1 processes. They do not 
contribute to the out of equilibrium condition. However they can lead to a dilution of 
the generated Ctotal- Therefore they are relevant in the exact computation of the n factor 
defined in Eq. (5.1). 

In order to compute the D — terms, we use the following relation between the equilib- 
rium densities: 



' ^fi{l-fl'){l-fn'){lTYc), 



(4.19) 



/Z(t)/f (1 + /|) = + /Z(t))(l + fn')^^'-'^ - 



(4.20) 



with 



feq _ \ 



(4.21) 



One gets 



Di + Di 



+ 



■fe|2 



•'jv 1 + ^ •'^ ■''^ ' 

■'Ni 



fk\2 



k 



where in order to write the equations in the closest to the standard notation we have 
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defined the following average widths: 

4' = - - /DElA'l' (4-22) 

k 

^iv,(4.) = ^i/iv.(l - - /DEl^l' (4-23) 



n 



^ "'at/ k 
n|(r|^.,) = A|/|(l + rJ){l + /^^) J^IA'I^ (4.25) 

^n.^n) = ^i/iv,(l + + (4-26) 



.eg/ps \ _ Al2 ^eg + ■^iV, ) 

which verify that in equilibrium 



4n,) = ^I/J'tttt^^i + -^F^^i + fH')E\^T\' (4-27) 

•'at/ k 



(r^f ) = (rgi) = (r^f ) (4-28) 

Equivalently for the rest of terms : 

A + A = - i;^(r^;) + ^:>^z:^/,(f (4.29) 

A - A = + )) + yi^y,(R ) (4.30) 

A + Dj = (y^^ - Yt^f^-}) + Yty^es, {%) (4.31) 
A - Dj = (y^^ (r|^^) + y|^(f + Y^JY,r~) (4.32) 

Concerning S" — terms, in order to evaluate, for example, the on-shell contribution 
\Mos{Lkh — Lkih)\'^ we use the following relation between the equilibrium densities 

(1 - )(i - /D = f~y^^^ [(1 - )(i - /D - /r] , (4.33) 

and the identity, 

1 = J d^N^^^^ -PL- Ph) ■ (4.34) 

They allow us to write the on-shell contribution to the scattering terms at the required 
order in e as: 

A?t /L/r(l - /l)(1 - fl')Y}HLkh ^ N,)\'^-—^Y.\^{N, L,rh)f = 



k "^Ni^ Ni k' 



X 



{2'Kf5'^{pj;r. -PL-Ph) X 



(27r)4 2m^r*^ 7 (27r)32£;L (27r)32£;ft ^ ' ^"^^^ 
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Now we define the thermal width into fermion and scalars as: 



k 

^iv, - 2mc.. 7 (27r)32£;.(27r)32EH^ ^ ^^^^ ^""^ 



(1 + + /h') - El^^l' • (4-36) 

fc 

so that the thermal width of sneutrinos is T*^ = T*-'-^ + r*~'^ 

Ni Ni Ni 

Using that / -p^27r5(p^^ - mj^) = J (_2ny^2E~ <=^" (^-3^) '^^s- 



2 r~'^ n ^2 r-'-'' 

N'r''' pt/i 



thj 



Altogether we find that 

S = Y^HLkh ^ Urh) - aiUrh ^ Lkh)) + y£e/,(r^^,,)-^ . (4.37) 

k,k' ' ' ATi 

The rest of on-shell contributions can be evaluated similarly: 

j-^th,s 

S = Y.{a{L,H lIh^) - a{LlH^ L,H)) + y|^e,,(r^.,)^ (4.38) 

k,k' ' ' Ni 



Sr7 



Ni 



k,k' ' JVj 

= ^{aiLkh ^ Lt,i/t) _ ^(Zt^i^t ^ i,/,)) + y£^!4±^(r^^,)_^ (4.39) 

= Y.{a{Uh ^ Lk,H) - a{Lk,H IJi)) - y^^!A±££i(r^^^)_^ (4.40) 

k,k' ' ' Ni 



' Ni 2 ^ Ar|=''/ rtfe 
k,k' Ni 

5^(a(Lfe/i ^ Lk,H) - a{L^,H L,h)) - YifJ^{r~^,)^ (4.41) 

k,k' ' ' Ni 
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fc,fc' 



Ni 2 



■IT'- 

\ AT 



N't"! Yth 



(4.42) 



Altogether we can write the Boltzmann equations for the sneutrinos and leptonic num- 
bers 



dt 



(4.43) 



dYc 
dt 



€f. ( Y^, (vt ) + Y^Vrt ) - 2Y^nrt\] + Y^'^Ycirt ) 

Tr \ Ni\ Ni' Ni^ Ni' Ni^ Nf' Ni Ni' 



-{2a{Lh Lh) + (7{Lh L^H^) + a{Lh LH)) 
+{2a{Lh Lh) + a{Lh LH) + a{Lh L^H^)) 
-{a{LH Lh) - a{LH Lh)) + 
+((7(Zti7^ ^ Lh) - u(pH^ Lh)) 



dt 



E 



[e,, { Yf-, {T% ) + YtVV- ) - 2YtVT%,,)) + Yt'^Y^iV- ) 

[ V N,' AT, \ Ni' AT, \ TVf'V AT, Ni' 



-{2a{LH L^H^) + a{LH Lh) + a{LH Lh)) 
+{2a{L^H^ LH) + a{L^H^ Lh) + a{L^H^ Lh)) 
-{a{Lh L^H^) - a{Lh LH)) 
+{a{Lh LH) - a{Lh L^H^)) 



(4.44) 



(4.45) 



The out of equihbrium condition is verified since using Eq. (4.28) the first term of Eq. (4.43) 

and the e terms of Eqs.(4.44) and (4.45) cancel out in thermal equilibrium. 

The Boltzmann equation for the total lepton number can be written as (here we use 

= -e/i = ei)- 



dt 



+ Y^{YcTt +Y-^T^~ ) + scattering terms 



+ s. t. 



(4.46) 



with (r^e,) = (r^,,) + (r^-^ 



^Here we have suppressed flavour indices in the two body AL — 2 scattering terms for the sake of 
simphcity, but a sum over afl flavours in initial and final states should be understood. 
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In the last line we have used that at 0{e) we can neglect the difference between 
and f- in the definitions of the thermal average widths and we have defined the effective 
T dependent total asymmetry: 

which in the approximate decay at rest takes the form Ref. [10, 11] 

e,(T) = 6,^^. (4.48) 

Cs + C/ 

5. Results 

We now quantify the conditions on the parameters which can be responsible for a successful 
leptogenesis. 

The final amount of — £ asymmetry Yg-/; = uq-c/s generated by the decay of 
the four light singlet sneutrino states Ni assuming no pre-existing asymmetry and thermal 
initial sneutrino densities can be parameterized as 

Yb-c = -i^Y. ^i^^d) Y'^XT » Mi) . (5.1) 

i 

ei{T) is given in Eq.(3.12) and is the temperature at the time of decay defined by the 
condition that the decay width is equal to the expansion rate of the universe: F = H{Td), 

1/2 T'S 1Q 

where the Hubble parameter H = 1.66 o* — — , m„; = 1.22 • 10 GeV is the Planck 
mass and ^r* counts the effective number of spin-degrees of freedom in thermal equilibrium, 
= 228.75 in the MSSM. Furthermore y!''(T > Mi) = 90C(3)/(47r''5*). 

In Eq. (5.1) k < 1 is a dilution factor which takes into account the possible inefficiency 
in the production of the singlet sneutrinos, the erasure of the generated asymmetry by L- 
violating scattering processes and the temperature dependence of the CP asymmetry ei{T). 
The precise value of k can only be obtained from numerical solution of the Boltzmann 
equations. Moreover, in general, the result depends on how the lepton asymmetry is 
distributed in the three lepton flavours [5]. For simplicity we will ignore flavour issues. 
Furthermore, in what follows wc will use an approximate constant value k = 0.2. 

After conversion by the sphaleron transitions, the final baryon asymmetry is related 
to the B — C asymmetry by 

UB _ 24 + Ann ub-c 
s 66 -|- I'inH s ' 
where nn is the number of Higgs doublets. For the MSSM: 



(5.2) 



-f = -8.4x10-^ Kj2^iiTd) (5.3) 

i 

This has to be compared with the WMAP measurements that in the ACDM model 
imply [27]: 

^ = (8.7;°i) X 10-11 (5.4) 
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Altogether we find that (for maximal CP violating phase sin(/> = 1): 

> 2.6 X 10-'^ (5.5) 



4 \BsN A\ r Cs{Td) - Cf{Td) ^ ^ ^ -,^-7 



4|S5ivP + |M|2r2 Cs{Ta) + Cf{Ta) 

Further constraints arise from the timing of the decay. First, successful leptogenesis 
requires the singlet sneutrinos to decay out of equilibrium: its decay width must be smaller 
than the expansion rate of the Universe T < H \t=m , with F given in Eq. (2.16), 



< 1.66 g^"— . (5.6) 



Stt nipi 
This condition gives an upper bound: 

Ein.i^(T)<^-»-' 

Second, in order for the generated lepton asymmetry to be converted into a baryon 
asymmetry via the B-L violating sphaleron processes, the singlet sneutrino decay should 
occur before the electroweak phase transition 



r > H{T ~ 100 GeV) ^ > 2.6 x 10"^^ GeV 

k 



(5.8) 



The combination of Eqs.(5.7) and (5.8) determines a range for the possible values of 
Yl l^ifeP foi" ^ given M: 

--0-(i5!^)<Em.l^<-10-.(^) ,5.) 

We now turn to the consequences that these constraints may have for the neutrino 
mass predictions in this scenario. Without loss of generality one can work in the basis in 
which Mij is diagonal. In that basis the light neutrino masses, Eq. (2.3), are: 

o s A V \2^^ lO^GeVlO^GcV 



where v = {H) is the Higgs vev. 

It is clear from Eq. (5.10) that the out of equilibrium condition, Eq. (5.7), implies that 
the contribution of the lightest pseudo-Dirac singlet ncvitrino generation to the neutrino 
mass is negligible. Consequently, to reproduce the observed mass differences AuIq and 
Am^j^, the dominant contribution to the neutrino masses must arise from the exchange 
of the heavier singlet neutrino states. 

This can be easily achieved, for example, in the single right-handed neutrino dominance 
mechanism (SRHND) [28]. These models naturally explain the strong hierarchy in the 
masses and the large mixing angle present in the light neutrino sector. In particular, in the 
simple case in which the matrix /j, and M are simultaneously diagonalizable, the results 
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Figure I: 2_^\^'i-k\'^ ~ s n regions in wliicli enough CP asymmetry can be generated (Eq. (5.5)) and 

k 

the non-equihbrium condition in the sneutrino decay, Eq. (5.6) and decay before the electroweak 
phase transition, Eq. (5.8) are verified. We take A = msusY = 10^ GeV. The regions correspond 
to M = 10^, 5 X 10"^, and 10^ GeV, from left to right. 



in Ref. [28] imply that for the inverse see-saw model with three generations of singlet 
neutrinos, the SRHND condition is attained if there is a strong hierarchy: 

Generically this means that the out of equilibrium condition requires the neutrino mass 
spectrum to be strongly hierarchical, mi <C m2 < m^. 

Conversely this implies that the measured neutrino masses do not impose any con- 
straint on the combination of Yukawa couplings and sneutrino masses which is relevant for 
the generation of the lepton asymmetry which can be taken as an independent parameter 
in the evaluation of the asymmetry. 

Finally we plot in Fig. I the range of parameters ^ l^ifcP and Bg^ for which enough 
asymmetry is generated, Eq. (5.5), and the out of equilibrium and pre-electroweak phase 
transition decay conditions, Eq. (5.9) are verified. We show the ranges for three values of 
M and for the characteristic value oi A = msusY = 10^ GeV. 

From the figure we see that this mechanism works for relatively small values of M 
(< 10^ GeV). The smaller is M, the smaller are the yukawas X^l^ifcP- Also, in total 
analogy with the standard seesaw [10], [11], the value of the soft supersymmetry-breaking 
bilinear Bsn, is well bellow the expected value MmsusY ■ The reason is that, in order 
to generate an asymmetry large enough Bsn ~ MF, but F is very small if the sneutrinos 
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decay out of equilibrium, F ^ 1 GeV (x(p^^^cv) ' ^ small B term with large CP phase can 
be realized naturally for example within the framework of gauge mediated supersymmetry 
breaking [29] or in warped extra dimensions [30] . 

Given the small required values of Bsn one can question the expansion in the small 
parameters in Eq. (2.13). As described at the end of Sec. 3 in order to verify the stability 
of the results we have redone the computation of the CP asymmetry keeping all the entries 
in the sneutrino mass matrix, and just assuming that it is real. We have found that as long 
as l-Bs'Arl 3> \Bs\, the total CP asymmetry is always proportional to Bsn, and presents 
the same resonant behaviour, so that it is still significant only for Bsn ^ MmsusY- 

In summary in this work we have studied the conditions for successful soft leptogenesis 
in the context of the supersymmetric inverse seesaw mechanism. In this model the lepton 
sector is extended with two electroweak singlet superfields to which opposite lepton number 
can be assigned. This scheme is characterized by a small lepton number violating Majorana 
mass term jj, with the effective light neutrino mass being nii, oc fi. The scalar sector contains 
four single sneutrino states per generation and, after supersymmetry breaking, their inter- 
action lagrangian contains both L-conserving and L- violating soft supersymmetry-breaking 
bilinear B-terms which together with the f^ parameter give a small mass splitting between 
the four singlet sneutrino states of a single generation. In combination with the trilinear 
soft supersymmetry breaking terms they also provide new CP violating phases needed to 
generate a lepton asymmetry in the singlet sneutrino decays. 

We have computed the relevant lepton asymmetry and we find in that, as long as the 
L-conserving S-term, Bsn, is not much smaller than the L-violating couplings, the asym- 
metry is proportional to Bsn and it is not suppressed by any L-violating parameter. As in 
the standard see-saw case, the asymmetry displays a resonance behaviour with the max- 
imum value of the asymmetry being obtained when the largest mass splitting, 2Bsn/M, 
is of the order of the singlet sneutrinos decay width, F. Consequently we find that this 
mechanism can lead to successful leptogenesis only for relatively small values of Bsn- The 
right-handed neutrino masses are low enough to elude the gravitino problem. Also, the out 
of equilibrium decay condition implies that the Yukawa couplings involving the lightest of 
the right-handed neutrinos are constrained to be very small which, for the naturally small 
values of the L-violating parameter n, implies that the neutrino mass spectrum has to be 
strongly hierarchical. 
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A. CP Asymmetry in Quantum Mechanics 



The four sneutrino system is completely analogous to the K° — K° system, so here we 
compute the CP asymmetry generated in their decay using the same formalism. In order 
to compare with the effective field theory approach described in Sec. 3 we consider only 
the simplified case B^, m^, m^, m|.jy <ti Bsn and ~ jiM* . In this limit, wc have 

chosen for simplicity a basis where A = |A|e"^ is the only complex parameter with (p given 
in Eq.(2.12). 

The evolution of the system is then determined by the effective Hamiltonian, 



H = M-i- 



(A.l) 



where, in the interaction basis and at leading order in the expansion parameters e, e defined 
in Eq. (2.13), 

/l e e\ 

1 e e 



M = M 



e e 1 
\i € 1 J 



(A.2) 



and 



/ 1 ^\ 



r = r 



1^0 
A 1 

^001/ 



(A.3) 



with r given in Eq.(2.16). 

It is convenient to write the effective Hamiltonian in the mass eigenstate basis Eq.(2.15), 
because in such basis the four sneutrino system decouples in two subsystems of two sneu- 
trinos, with the resulting width matrix: 



r 



f 1 — eA cos ( 
eA sin (f) 


V 



CA sm (p 
1 + CA cos (f) 
1 







eA cos ( 






-eA sin( 



-eAsincf) l + eACOS(pJ 



where = pj. 



The eigenvectors of the effective Hamiltonian H are: 



(A.4) 



2+2 



If+I 



2+2 



If+l 



2+2 
1 



2+2 



If+I 
l^-l 



gi(,^_-<^+)/4(^^t _iVt) + ^ 






-0+)/4(5 _ TV) 


gi(0--<^+)/4(5t_iVt) - . 


fHe- 




-'^+)/4(5_ JV) 


gi(<^_-<^+)/4(^t+ivt)_ 


1 k-l^ 




-<!>+) (S + N) 


ei(<A--<^+)/4(5t +ivt) + . 


/ k-l^ 




-'t>+)/i(S + N) 
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and the eigenvalues 



1^4 



where 



: \M\ - 1^1/2 - ir/2 - e^('^-+'^+)/2y^y^ 
: \M\ - 1^1/2 - ir/2 + e^('^-+'^+)/2y^y^ 
: |M| + 1^1/2 - ir/2 - e^('^-+<^+)/2y^y^ 
: |M| + I//I/2 - ir/2 + e^('^-+^+)/2 v^^T^ 

e_ = |e_|e^<^" = |B57v/(2M)| - ireA/2 
e+ = |e+|e^<^+ = |Bs7v/(2M)| - iVe\l2 

We consider an initial state at t = with equal number densities of the four sneutrino 
interaction states F^. Using that the time evolution for the hamiltonian eigenstates Eq.(A.5) 



is trivially 
following: 



,—ivit 



N-), we obtain that at time t the interaction states are the 



N{t) 



gi+it) + g2+it) 



N) + e^('^— <^+)/2 9i-{t) -92-{t) 



9i+{t) - g2+{t) 



S) - e'{4>--4>+)/2 gi-{t) +92-{t) 



5t 



e+[ i(^__</,+)/2 9i-{t) - 92-{t) 

4 



N) + 



gi+it) + g2+{t) 



/ k+l -ii^—<p+)/2 gi-{t) + g2-{t) 

\ le_| 4 



s)- 



gi+it) - 92+{t) 



5t 



N 



iVt\ + 



+ 



5i+(t) + 52+(i) 



5t 



-i(^--<p+)/2 gl-W +ff2-(t) 

Ie_| 4 



TV - 



gi+{t) - g2+{t) 



iVt 



+ 



_-0+)/2 gi-{t) - g2-{t) 

4 



5 + 



ffl+W+52+(t) 



5t 



The functions gi±{t) and 52±(i) containing the time dependence are given by: 



(A.5) 



(A.6) 
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with, 



(A.7) 



(A.8) 



We neglect soft supersymmetry-breaking terms, so that S only decay to scalars, and 
N to antifermions: 

^ \A[S L,Hr = \A[S^ ^ lIh^I' = \Y,,M\' ^ \A~^f 
\A[N Llh^\^ = \A[N^ ^ L,hr = \Yr,\Hs - ml - ml) ^ \ALf ^ ' ' 

We next write the time dependent decay amplitudes in terms of I^J^P p: 



\A[N{t)^L,H]\' = \A[N\t)^LlH^\ 



to-tl|2 _ - 



16 



(A.IO) 



\A[N{t) LlH^f = (M) V[Art(t) ^ = j^^ \A~^f{AAl) 

\A[N{t) L,h]f = V[iVt(t) ^ 4/,t] = |!^ \9i--92-\' |^^j2(A.i2) 



't J,tl|2 _ I /lr^rt^^^ , r_f,ll2_ lfl++52+P 



tf7-tl|2 _ 



\A[S{t) LlH^]f = (1^) \A[S\t) LkH]f = |^ ^/^"^ lA^p (A.15) 
\A[S{t) ^ = (^My|^[5t(t) ^ 4/jt]|2 = 1^^-+/^-!' \AlJ\AA6) 



\A[S{t)^Llh^\' = \A[S\t)^L,h]\' = 



16 

51+ + 52+1 ' 



I^L.P (A.13) 



16 



1^ 



2 



(A.14) 



51+ -52+1 ' 



16 



(A.17) 



We define the integrated CP asymmetries for the fermionic and scalar channels as: 

/ dtj2 [\A[m) -^Lk+xr- \A[m) ^ 4 + xjpi 



i,k 



i,k 

jdtY^ [\A[m - L,+xr - \A[m ^ ii+xf 

i,k 



jdtY^ [\Am) ^ + + \Amt) - 4 + xf 

i,k 



(A.18) 



Using the time-dependent amplitudes from Eq. (A), the time integrated asymmetries 



are: 



-e/ = e 



1 / Ie 



k+l 

2Vk+| 



(A.19) 
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where the factor |^ — vanishes if CP is conserved, and for ^ 1 is given by: 



|e_| |e+| 4Mr|eA| 4ryl , 

M-^- „ "^' sin<^=-— sin<^ (A.20) 
|e+| |e_| BsN osN 



The time dependence is encoded in x'- 



^ /~di[|5i+(i)P + l52+(i)P + l5i-WP + l52-(*)P] ■ ^' ^ 



In the hmit ca ^ 1 the time integrals are : 

poo o 

dt + b2+p + + 152- (Oil =^ r ' ^^-^^^ 

The comparison between the asymmetry in Eq.(3.11) and Eq.(A.24) is in full analogy 
to the corresponding comparison in the standard see-saw case discussed in [10] . The asym- 
metry computed in the quantum mechanics approach, based on an effective (non hermitic) 
Hamiltonian, Eq.(A.24), agrees with the one obtained using a field-theoretical approach, 
Eq.(3.11) in the limit V <C Bsn/M. When T Bsn/M, the four sneutrino states become 
two pairs of not well-separated particles. In this case the result for the asymmetry can 
depend on how the initial state is prepared. If one assumes that the singlet sneutrinos 
axe in a thermal bath with a thermalization time F"^ shorter than the typical oscillation 
times, AM^^, coherence is lost and it is appropriate to compute the CP asymmetry in 
terms of the mass eigenstates Eq.(2.15) as done in Sec. 3 and one obtains Eq. (3.11). If, 
on the contrary, one assumed that the N, S states are produced in interaction eigenstates 
and coherence is not lost in their evolution, then it is appropriate to compute the CP 
asymmetry in terms of the interaction eigenstates Eq.(A.5) as done in this appendix and 
one obtains Eq.(A.24). 
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